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Abstract: We consider a GI/G/c/K- type retrial queueing system with constant retrial 
rate. The system consists of a primary queue and an orbit queue. The primary queue has c 
identical servers and can accommodate the maximal number of K jobs. If a newly arriving 
job finds the full primary queue, it joins the orbit. The original primary jobs arrive to the 
system according to a renewal process. The jobs have general i.i.d. service times. A job in 
front of the orbit queue retries to enter the primary queue after an exponentially distributed 
time independent of the orbit queue length. Telephone exchange systems, Medium Access 
Protocols and short TCP transfers are just some applications of the proposed queueing sys- 
tem. For this system we establish minimal sufficient stability conditions. Our model is very 
general. In addition, to the known particular cases (e.g., M/G/l/1 or M/M/c/c systems), 
the proposed model covers as particular cases the deterministic service model and the Erlang 
model with constant retrial rate. The latter particular cases have not been considered in 
the past. The obtained stability conditions have clear probabilistic interpretation. 
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Stochastic Stability 
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L'analyse de la stabilite de la file d'attente de type 
GI/G/c/K avec des clients qui reviennent a un taux 

constant 

Resume : On considere une file d'attente de type GI/G/c/K avec des clients qui re- 
viennent a un taux constant. Le systeme se compose d'une file d'attente primaire et une 
file d'attente orbite. La file d'attente primaire a c serveurs identiques et peut accueillir le 
nombre maximal de K clients. Si un arrive trouve la file d'attente primaire pleine, il re- 
joint l'orbite. Les clients qui entrent dans le systeme pour la premiere fois arrivent selon un 
processus de renouvellement. Les clients ont un temps de service generale iid. Les clients 
dans la file d'attente orbite essaient d'entrer dans la file d'attente primaire apres un temps 
avec une distribution exponentielle independante de la longueur de la file d'attente orbite. 
Les commutateurs telephoniques, le controle d'acces au support, et les courte transferts TCP 
sont quelques-unes des applications de le systeme etudie. Pour ce systeme, nous etablissons 
les conditions de stabilite suffisantes. Notre modele est tres general. En plus des cas parti- 
culiers (par exemple, M/G/l/1 ou M/M/c/c) , le modele propose couvre les cas particuliers 
du modele de service deterministe et le modele Erlang avec des clients qui reviennent. Les 
derniers cas particuliers n'ont pas ete considere dans le passe. Les conditions de stabilite 
obtenus ont une interpretation probabiliste tres claire. 

Mots-cles : File d'Attente avec des Clients qui Reviennent, Stabilite Stochastique, Proces- 
sus de Renouvellement, Regeneration 
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1 Introduction 

We consider a GI/G/c/K-type retrial queueing system. The system consists of a primary 
queue and an orbit queue. The primary queue has c identical servers and can accommodate 
the maximal number of K jobs. If a newly arriving job finds the full system, it joins the orbit. 
The original primary jobs arrive to the system according to a renewal process with rate A. 
We denote the arrival times of the original primary jobs by t n and we denote the interarrival 
times by r„ = £ n +i — t n , n > 1, with generic element r. Without loss of generality we assume 
that t\ = 0. The jobs have general i.i.d. service times {Sn , n > 1} with service rate \i 
and generic element S^ 1 ' . Retrial times {S n , n > 1} are i.i.d. exponential with (generic) 
service time 5*^ 2 ^ with rate /j,q and independent of the orbit size, provided it is positive. Such 
a model is referred to as a retrial model with constant retrial rate. It then follows that the 
orbit can be interpreted as a single-server -/M/l-type queue with service rate /xo and with 
the input which is formed by the flow of jobs rejected from the primary queue. We use 
notation -/M/l because the merged stream arriving in the orbit is not in general G/-type 
since it is a complex combination of a lost part of original primary customers and secondary 
customers returning to the orbit after unsuccessful attempts to enter the primary queue. 

In this work we establish minimal sufficient stability conditions for the presented retrial 
system. Our model is quite general. To the best of our knowledge, the retrial queueing 
system with constant retrial rate and general renewal arrival process is considered for the 
first time. In [TT] Fayolle has introduced a retrial system with constant retrial rate. Fayolle 
has derived stability conditions for the case of M/G/l/1 primary queue. In Q] Artalejo 
has obtained stability conditions for the Markovian case M/M/2/2. In [20] Ramalhoto and 
Gomez-Corral have obtained stability conditions for the M/M/l/2 case. For the general 
Markovian case M/M/c/K the authors of [50] have obtained decomposition results assuming 
ergodicity. The ergodicity conditions for the multiserver Markovian case M/M/c/c with 
recovery probability have been derived by Artalejo, Gomez-Corral and Neuts in [2J. In [T3] 
Krishna Kumar and Raja have derived stability conditions for the M/M/c/c constant retrial 
rate model with feedback and balking. Stability conditions for the basic M/M/c/c constant 
retrial model can be recovered from the results in [5] and [T3] . The ergodicity conditions for 
the Markovian case M/M/l/K have been obtained in [2]. In the works [2J, [I] and P2] the 
authors have established stability conditions by using the matrix-analytic technique for QBD 
processes [TH]. We shall demonstrate that stability conditions for M/G/l/1, M/M/c/c and 
M/M/l/K systems are particular cases of our general conditions. Furthermore, the following 
important cases with Poisson input have not previously been covered: M/G/c/c (Erlang 
model), M/M/c/K (Markovian model with a general number of servers and waiting spaces), 
M/D/l/K system. Stability conditions of these important cases appear to be particular 
cases of our general conditions and will be considered in detail below. We emphasize that, 
to the best of our knowledge, such retrial system with general renewal input of primary 
customers is considered for the first time. 

There is a number of applications of retrial systems with constant retrial rates in telecom- 
munication. Using a retrial queue with constant retrial rate Fayolle |11| has modelled a 
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telephone exchange system. In the series of papers 0-111!] the authors have proposed to use 
a retrial queueing system with constant retrial rate to model Multiple Access protocols. In 
particular, in [S] the authors have modelled an unslotted Carrier Sense Multiple Access with 
Collision Detection (CSMA/CD) protocol and in [S] and [TO] the authors have modelled 
some versions of the ALOHA protocol. In [JJ and [5] the authors have suggested to use re- 
trial queues and retrial networks with constant retrial rates to model TCP traffic originated 
from short HTTP connections. 

The stability analysis used in this paper is based on renewal theory and a characterization 
of the limiting behavior of the forward renewal time in the process generated by regenerations 
of a basic process. This approach, presented in a general form in [161 118) . turns out to be 
effective in the stability analysis of many queueing systems including general multiserver 
retrial queue |17| , and also multiserver system with non-identical servers |15| . The presented 
method also works successfully outside of Markovian models, and it is demonstrated also in 
this paper where non-Markov processes are considered. In particular, it allows us to reduce 
the dimension of the processes that simplifies analysis considerably and does not require 
involved stability techniques developed in the theory of multi-dimensional Markov processes 
[5J I14| . An important contribution of this work is an extension of the stability analysis to 
arbitrary initial state of the system. 

The paper is organized as follows. In the next Section 2 we present the main result of the 
paper, minimal stability conditions for GI/G/c/K-type retrial queue with constant retrial 
rate. We also provide the proof of the main result. Then, in Section 3 we specify the general 
condition for a number of important particular cases of the primary queue. In particular, we 
consider the general Markovian queue, the Erlang queue and the queue with deterministic 
service. We would like to note that stability conditions for these important particular retrial 
queues have not been available before. We conclude our paper with Section 4. 

2 Stability analysis 

To describe the behavior of the system, we consider (right-continuous) process M(t) :— 
N(t) + v{t), t > 0, where N(i) is the number of retrial customers being in orbit and v(t) is 
the number of the customers waiting in the buffer, at instant t. Note that v(t) £ [0, K — c] 
for any t. Also we introduce the (right-continuous) process W(t), t > 0, expressing, at each 
instant t, the remaining workload in all servers. More exactly, if Si(t) is the remaining 
service time at server i at instant t (= if the server is empty), then W(t) = Xa=i Si(t). 
Introduce the basic (two-dimensional) process X = {X(t) := (M(t), W(t)),t > 0}. The 
choice of the basic process is motivated by the stability analysis under arbitrary initial state 
of the system, while there are other candidates for basic process, to analyze the system 
under zero initial state. In the latter case, one can use, for instance, the total number of 
customers in the system, or the total remaining workload in the system (including orbit), 
etc. Also, let X(t^) = Xk, k > 1. Denote To = 0, then the instants 

T n+1 = inf (t k > T n : X k = (0, 0)), n > 0, (1) 
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are the regeneration points of the basic process X. Let T be a generic regeneration period 
and T(t) = inffc{Tfc — t : Tk — t > 0} be the forward renewal/regeneration time at instant 
t > 0. If 

T\ < oo with probability 1 (w.p.l) and ET < oo, (2) 

then we call any regenerative process having regeneration instants defined by (TTJ) (and also 
the original system) positive recurrent. (This term, repeatedly used in previous works [151 
[TCI [TT1 \Ti\ \TI\ . has evident analogy with positive recurrent Harris Markov chains whose 
embedded renewal process of regenerations has finite mean cycle length.) 

It follows from the theory of regenerative processes that positive recurrence is the most 
essential element in stability analysis of the process. Indeed, if the interarrival time r is non- 
lattice then T is so, and, under positive recurrence, X{t) converges to a stationary limit. (It 
is obvious that the system is unstable if ET = oo.) Our approach to stability is based on 
the following result [12]: if ET = oo then 

T(t) => oo as t ->■ oo, (3) 

regardless of initial (finite) delay T\ = T(0), where =>■ stands for convergence in probability. 
Thus, if convergence ((3|) does not hold then (0) is satisfied and positive recurrence for the 
zero- delayed process holds, in which case T\ = T and X(Q) = X\ = (0, 0) (that is t\ = is a 
regeneration instant). However, violation of ([3]) does not implies in general the finiteness of 
T\ w.p.l under arbitrary initial state X(0), and an extra (sometimes hard) work is required 
to prove that T\ < oo. Denote fto = 0, then 

0n +1 = in£(k>p n :X k = (0,0)),n>0, (4) 

k 

are regeneration instants of the embedded discrete-time process X n , n > 1, with generic 
regeneration cycle length /3, and its remaining renewal time j3(n) := iufk{(3k — n : f3k — n > 
0} => oo (as n — » oo) provided E/3 = oo. 

Now we outline how to apply the above mentioned approach to our model. Because the 
buffer of the primary queue is finite, the source of instability of the system can only be 
the unlimited increase of the orbit size. Thus, we first show that under the predetermined 
stability condition (see Theorem 2.1) the orbit size N(t) ^4> oo. Then the second step is to 
show that also T(t) =4> oo. Finally, we apply characterization ©■ To establish T\ < oo for 
non-zero initial state, we will use new development of the approach presented in |18| . 

Let us denote original system by £ and construct an auxiliary (new) system S as follows. 
The system E has the same set of servers and the same buffer as system E, the same renewal 
input of primary customers with rate A (we call them A-customers) and in addition an 
independent Poisson input of primary customers with rate [Iq (we call them /zo-customers). 
Arriving primary customer (of any type) who finds servers and the buffer full in system 
E joins the orbit. The orbit is an infinite buffer system of ■ jMj\ type as in the original 
system. The secondary customers leaving orbit (in E) leave the system forever and do not 
affect its future state. Note that the system E (as E) regenerates at the instants when the 
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A-customers find both buffer and orbit empty. Note that here we also use the memoryless 
property of the input of /^-customers. Moreover, for any variable C, in system E we denote 
corresponding variable in system E as £. In particular, W k is the remaining workload in all 
servers and v k is the number of waiting customer in the buffer, respectively, in the system 
£ at instant tZ, k > 1. 

Denote the primary queue of system £ as 0, and note that it can be considered as an 
isolated system because secondary customers (in £) leave the system and do not go to the 
primary queue. Then the subsystem regenerates (in continuous time) at the instants 

*„+i = mm{t k > #„ : Wj. = z> fe = 0), n> 0, (5) 

where, by definition, \&o = 0. (So the regenerations of the whole system E are a subsequence 
of the regenerations of 0.) Because the subsystem has a regenerative input and finite 
buffer, then such a system is positive recurrent under the condition 

P(r>5 (1) )>0, (6) 

see [161 118] . In particular, mean generic regeneration period E\I/ < oo. 

Remark 1 Exact condition in J 16!/ (condition (3.15) there) and in ]!$/ (condition (31) 
there) require that the discrete-time regeneration period A (counting all arrivals during re- 
generation period $>) equals 1 with a positive probability. It occurs if the interarrival time 
t following a X-customer (with service time S^ 1 ' ) starting new regeneration period is larger 
than S^ 1 ' and less than the next Poisson arrival (with interarrival time f). In other words, 
the following inclusion holds: {A = 1} D {f > f > S^ 1 '}. Obviously, under condition $j\), 
P(f > t > S^ 1 ') > 0, and the required assumption is fulfilled. 

Denote by R(t) the total number of rejected customers in system in the interval [0, t] 
(this is also the total number of customers which went to the orbit in the whole system E). 
Denote by A(t) the total number of arrivals (primary A-customers and /xo-customers) in the 
interval [0, i\. Denote also by R the number of rejected customers during regeneration cycle 
of system 0. 

Of course, the process {R(t), t > 0} is positive recurrent cumulative process with em- 
bedded regenerations {^n} and, in particular, there exists (w.p.l) the limit 

hm ^1 = IS. (7) 

*^~ A{t) EA 

(To explain, we note that R(t)/t ->■ E^/E*, A(t)/t -)■ E.4/E*.) In the system E, define 
indicator /„ as 

1, if customer n is rejected, . . 

0, otherwise, ^ ' 

so the sequence {/„, n > 1} has regeneration period A. Because P{A = 1) > 0, then the weak 
limit I n => / exists, or P(/„ = 1) — > EI := Pi OS s where Pi oss is stationary loss probability. 
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Moreover, by the standard result of regenerative theory, stationary loss probability coincides 
with the long-run-average loss probability (Q : 

P loss = ton ^ = ^. (9) 

Now we are ready to formulate the main stability result. 
Theorem 1 Assume that condition © and the following condition 

{\ + l-lo)Ploss < Mo, (10) 

hold. Then, under arbitrary (fixed) initial state X(0) = X± = (mo, Wo), the original system 
is positive recurrent, that is 

ET < oo, E/3 < oo and T\ < oo, /?i < oo w. p. 1. 

Moreover, the stationary distribution lim^-nx; P(X n £ •) exists. If, in addition, interarrival 
time t is non-lattice, then stationary distribution ton.t_j.oo P(X(t) G •) also exists. 

Proof. We will use a monotonicity property of the loss system with respect to change of 
service times, see [2H [551 [__] ■ Namely, we use coupling to sample identical corresponding 
interarrival times in both systems £, £ until first empty period of the orbit in the original 
system £. Then we assume that during this (and any following) empty period of the orbit 
we continue to sample input Poisson process (with parameter no) in the original system 
but with arriving customers having zero service times (as long as the empty period lasts). 
At the end of such period, we interrupt current interarrival times (in Poisson processes in 
both systems £ and E) and resample new (identical) interarrival times for both systems. 
By memoryless property, this resampling keeps distribution of the input process. At the 
same time this procedure allows us to keep equivalence between input intervals in both 
systems. (Note that A-inputs stay identical and unchanged in both systems.) Moreover, this 
shows that service times of arriving customers in the systems are stochastically ordered as 
Sn <st Sn while for non-zero (actual) customer n, Sii ~ st Sn , n > 1. 

Now we can use the result of [23] which claims (in adaptation to our model) that if 
two (the same) finite capacity systems have the same input and the ordered service times, 
as above, then the number of rejected customers in primary queues of systems E and E, 
respectively, in the interval [0, t] are ordered as 

R(t)< st R(t), t>0. (11) 

We emphasize that R(t) counts all rejections which happen in system E including repetitive 
rejections of the orbit customers after unsuccessful attempts to enter the primary queue. 
Because the input in system E is a superposition of two independent renewal processes, 
then w.p. 1 as t — > oo 

^^X + fio- (12) 
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Since R(t) = X^ (t) 4, it then follows from © and © that w.p.l 

l im M = ( A + Mo )p, oss . (13) 

t— s-oo t 

(It follows from #(*) < A(t) and EA(i)/i -> A + fi , that the family {R(t)/t, t > 0} is 
uniformly integrable, and thus the limit ER(t)/t — >• (A + /io)Pz OS s also exists.) 

Denote by V (£), V^,(i) the total workloads arrived to the orbit in the systems S, £, 
respectively, during the interval [0, t], including the same (arbitrary) initial state V o (0) = 
V o (0) :=V . Note that 

R(t) 

V (t) = V o (0) + ]T Sj?\ t>0 (£ = 0), (14) 

fc=l 

and thus, 

tw i , . , 1K , 

; ->■ — w.p.l as i -)■ oo. (15) 

i?(i) Mo 

Denote, in the original system E, by /u (i) the total empty time of the orbit in interval [0, t] 
and by W (t) the (right-continuous) remaining workload in orbit at instant t > 0. Now we 
have the following balance equation 

Vo(t) = W (t)+t-iio(t), t>0. (16) 



Thus, we have 



Ho® > t - V (t) =t-J2 4 - VoQS) 



R(t) 

(2) 



fe=l 
R(t) 



> 



fit 



1 - Yl ^ - ^(°) = l - ^"(t), *>0. (17) 



fe=i 



Bv ([ID] ) . ([T5 ]) . ([15 ) this implies 

l iminf MoW > ! _ (A + ^)PlQ M ;= 5q > (18) 

This also shows that under arbitrary initial state V o (0), 

Mo(*) = / I(N(u) = 0)du — > oo w.p.l as i — > oo, (19) 

Jo 
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where / denotes indicator function. By Fatou's lemma, 

liminf-E/u (i) > 0. (20) 

t— >oo t 

Thus, P(N(t) = 0) /> as t — > oo and there exist 6* > and non-random sequence of 
instants z n — > oo such that 

inf P(N(z n ) = 0) > 5*. (21) 

n>l 

Starting from this point we work with the original system E only. Let r(t) = minfc(^ — t : 
tk ~ t > 0) be the remaining interarrival time at instant t in system E. Note that, for each 

P(N(z n ) = 0) = P(jV(z„ +r(z n )) = 0) > 5*. (22) 

Denote k(u) = min(fc : t^ > z n ), then z n + r(z n ) = t K („) is the first arrival instant after 
z„, n > 1. Let L(t) be the total remaining workload at the servers and buffer at instant t. 
Obviously, 

K-c 

m^tW^+^S^, t>0. (23) 

i=\ 

Note that the remaining service time processes, {Si(t), t > 0}, i = 1, . . . ,c, are tight. (To 
have the tightness, the only requirement, besides the ffniteness of EiS**- 1 ), is that service times 
{Sn } are independent of the input process [IS]-) Thus, the process W(t), t > 0, is tight, 
and hence, the process L(t), t > 0, is also tight. Then we can find a constant D such that 
inf„ P(L(z n ) < D) > 1 — 5* /2. Now we fix for a moment some z n satisfying (|2"Tj) and let 
•^(^n) = ^ni n. > 1. Then, we have 

P(7V K( „) = 0, L K{n) <D)> P(N(z n ) = 0, L(z n ) < D) > y, (24) 

where the first inequality is valid since the workload decreases in interval [z n , t K r n \) and the 
second inequality follows from 

P(N(z n ) = 0, L(z n ) <D) = P(L(z n ) <D)- P(N(z n ) > 0, L(z n ) < D) 

> P(L(z n ) <D)- P(N(z n ) > 0). 
Next we introduce the event 

£n ■= {N K{n) = 0, L K(n) < D, r(z n ) < C} 

where, by the tightness of the process {r(t),t > 0}, the constant C is taken in such a way 
that P(£ n ) > <5*/4. Thus, on the event £ ni the customer k(ji) indeed arrives in interval 
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[z n , z n + C] and finds the workload L K i n \ < D. It follows from ^ and Et < oo that one 
can find constants a < oo, 7 > 0, e > such that 

P(a>r >j + S {1) ) = e. (25) 

Note that, for each i, on the event w, := {a > n > 7 + S^ } the workload accumulated 
in the system at instant ti, decreases during interarrival time [£,-, £j+i) not less than by 7, 
provided that at least one server is not empty during this interval. If only orbit is not empty, 
then an orbit customer may attempt to enter server/buffer before next arrival instant U+i. 
In such a case, regardless of the retrial attempt was successful or not, on the event u>i, the 
accumulated workload decreases not less than by 7/2 (during [£$, i«+i)) with the probability 
> e(l - e^"2). Denote R := \2D/-/]. Then on the event 

i— K,(n) 

R primary customers arrive in the interval [z n , z n + C + aR], and the accumulated workload 
reaches zero (regeneration) with probability > j6*e R (l — e~^°^) R > 0. Because this lower 
bound is uniform in n, z n , then E/3 < 00. It follows from Wald's identity and representation 
T =st y^ft._i Tk that mean cycle period for continuous time processes is finite, ET = E/3Er < 
00. 

Using a modification of the approach from |18) . we now extend analysis to arbitrary 
initial state of the basic process X. Let S k ,K be the fc-th service time realized at server i. 
(For each i, these service times are i.i.d. and distributed as S^\) Introduce 

Si(t) = min^ + • ■ • + Sg } - t : S™ +■■■ + 5« - t > 0) 

the remaining renewal time at instant t in the renewal process generated by service times 
of server i, and let Si(0) = Si(0), i = 1, . . . , c. For any (integer) r > and x > 0, denote 
the set B(r, x) = [0, r] x [0, x] and consider the process Y =: {Y(t) = (N(t), W(t)), t > 0}. 
Then, for any t and x: 

I(Y(t) S B(0, x)) > I(N(t)=0)-I(W(t)>x) 

C 

(26) 



Although, for each t, Si(t) and Si(t) in general are not comparable, but by construction, 
the following inequality holds for any x: 



t , ft 



l(Si(u)>x)du< l(Si(u)>x)du, i = l,...,c. (27) 
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Note that l(§i(t) > x], t > 0, is a regenerative process with (generic) cycle length S^\ 
Then w.p.l as t — > oo, 

j f l(Si(t) > x) ^ ^ f°° P(S« > t)dt := P(SP >x), i = 1, • . . , c, (28) 



ii 



where distribution of the stationary overshoot Se is proper. Hence, we can take x\ in such 
a way that 

P(5« > ^) < |. (29) 

If we take x = x 1 in $Z$-$g% then it follows from (JTHJ), ©-flU) that w.p.l 



1 f* 

liminf- / I(Y(u) G B(0, Xi))du > liminf 



MqW _ ^o _ So ^ 



In particular, the total time that the process Y spends in the set B(0, x\) within time interval 

[0,t] 



/ I(Y(u) G B(0, xx))du ->• oo, t -> oo. 
Jo 



(31) 



Because, for each £ and x > 0, 

{Y(t) G B(0, x)} C {X(t) G B(K - c, x)}, (32) 

then, by (f3"Tj) . the total time the process X spends in the set B(K — c, xi) during time 
interval [0, t] 

HB(t) '■= I I{X{u) G B(K — c, xi))du -4 oo w.p.l as t — ¥ oo. 

JO 

Denote by A(t) the number of A-customers arriving in the interval [0, t], and let Gs(i) be 
the number of these A-customers who meets the process X in the set B(K — c, xi), that is 

A(t) 
G B (t) = Y, HXk G B(A' - c, xr)). (33) 

fe=i 

It then follows that w.p.l 

Gs(t)( max r„) > Ub(*) — > oo, t — > oo. (34) 

l<n<A(t) 

Since Er < oo, then maxi<„<A(() t„ = o(t), t — > oo [23] , Thus, Gs(£) — > oo w.p. 1 as 
£ — > oo, and the total number of A-customers, Gs(oo), meeting the process X in the set 
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B(K-c, xi) is infinite w.p. 1. (This holds for arbitrary initial state X(0), see ([Ig ]) .([Tg ]) .([3"0 |) .) 
Recall initial state X(0) — (mo, wq), and choose and fix arbitrary r > max(?7io, K — c) and 
z > max(wo, Xi). Then, as in |18j . one can find constants Dq < oo, £o > such that the 
remaining renewal time f3{k) satisfies 

inf P {mo , Wo )(P(k) <D \X k e B(r, z), fa > k) > e , (35) 

where the bounds Sq and Dq are uniform: they depend on r, z but does not depend on 
neither customer number k nor the specific states Xk e B(r, z). Indeed, due to the specific 
form of the basic process X, conditioned on the event {Xk £ B(r, z), Pi > k}, the total 
workload in the system at instant tZ is (uniformly in k) stochastically upper bounded by 
the quantity 

s+eV+e^, (36) 

where, evidently, both the service times \S\ } of waiting customers in the buffer and (expo- 
nential) retrial times \S\ } of all orbit customers are i.i.d and independent of the remaining 
workload Wk (< z) in the servers. Then one can unload the system during Dq A-arrivals 
with a probability > eo (like in the proof after formula ([3H|) ). and thus condition (pJ5"| is 
indeed satisfied. Then one can show that the mean number of A-customers meeting the 
process X in the set B(r, z) during first regeneration period is finite, 



t(m ,«>o) 



(E^fcGB(r, z)))<^<oo. (37) 



(More details can be found in [TTl [T5] .) Because B(r, z) D B(K — c, xi), then ([3T]) implies 
that the number of A-customers meeting the process X in the set B(K — c, x\) during first 
regeneration period is finite, G B := Y^k=i ^(^fc e B(/\ — c, x\)) < oo. Then it follows from 
above that P( mo ,w )(^'B ^ ^b(oo)) = 1, and hence the number of regeneration cycles is 
not less than two. Thus, first regeneration period pi < oo and hence T\ = T\ + • • •+Tp l < oo. 
(Indeed, the number of regeneration cycles is infinite w.p.l.) Finally, we note that condition 
(O implies aperiodicity of ft. ■ 



Remark 2 In general, one can define another type of regeneration, for instance, k-regeneration, 
when an arriving X-type customer see all servers and buffer empty and a fixed k customers 
in the orbit, see /J7f . Moreover, a relaxation (or elimination) of assumption (0|) (which 
is necessary to have classical regenerations of type {I]), H\)) may lead to the so-called one- 
dependent (or weak) regeneration with a dependence between adjacent regeneration cycles, 
see for instance, ]1 61 \21l \2^ . However, we do not consider such regenerations in this work. 
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Remark 3 In a quite general case of unbounded t (and in particular for the Poisson X- 
arrival process), the constant a in h2b}) can be taken in such a way that regeneration occurs 
in the interval [z n , z n + C + a]. 

3 Application to particular queueing models 

In this section we apply general stability conditions of Theorem [1] to important particular 
types of the primary queue. 

3.1 Markovian case M/M/c/K 

In the Markovian case with c servers, service rate /i and K places in the primary queueing 
system, the loss probability in the auxiliary £2 system is given by 



where 



Po = 



E 



((A + W))/^) c /(A + /i )//i 



K-c 

E 



K- 



Po, 



((A + mo)/m)" , ((A + ^o)//i) c 



(A + Mo)//" 



Hence, the stability condition (jTU|> takes the form 
((A + mo)/m)* ({X + noVn 



K-c 



Pn< 



/'o 



A + Mo 



(38) 



The other conditions of Theorem Q] are naturally satisfied by the Poissonian arrival process. 
In the case of one server (M/M/l/K case), the above condition reduces to 



((a + Mo )/m) j 



K 



£((A + M o)/M) ? 



< 



/'o 



A + Mo' 



which is equivalent to the ergodicity condition provided in 0]. Then, in the M/M/c/c case, 
the condition (|38|) reduces to 



((A + Mo )/m) c 



E 



((A + Mo )/M) r 



< 



/'o 



A + Mo' 



(39) 



which is equivalent to the ergodicity condition provided in [5] with the recovery probability 
p=l. 

It is important to note that the conditions given above are in fact necessary and sufficient 
stability conditions for the Markovian case. This is why we call conditions presented in 
Theorem [T] minimal sufficient stability conditions. 
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3.2 Erlang model 

When the primary queue is described by the Erlang model (M/G/c/c case), the general 
stability condition (jTUJ) takes the form 



((A + Mo )ES« 



(A + /x )EgW) 



-1 



< 



A + Mo' 



which is the same form as in (|55|) . Thus, the stability conditions demonstrate insensitivity 
property when the primary queue is the Erlang queue. We would like to emphasize that it 
seems to be very difficult (if possible at all) to obtain the above stability condition using 
stability techniques for Markov chains and embedded Markov chains. Anyway, the presented 
approach ensures short and simple proof of stability. 

3.3 Deterministic service 

We can obtain explicit stability conditions in another important case, in the case of de- 
terministic service (M/D/l/K model for the primary queue). In the case of deterministic 
service, the loss probability is given by 

P, -1- hK - X 

± loss *• 



1 + pbic-1 

where 

Y { -^-{K - 1 - n) n p n e^ K - 1 - n ^, 



K-\ 

!*-! = 



n=0 

with p = XS j7]. Thus, in the M/D/l/K case the stability condition (|TU|) takes the form 

A b K -i 



A + Mo 1 + pbx-i 
and we note that there exists an efficient recursive approach to calculate the coefficient bx-i 

4 Conclusion 

We have considered a retrial queueing system with general renewal arrival process, general 
service time and constant retrial rate. To the best of our knowledge a retrial queueing system 
with constant retrial rate was considered for the first time under so general setting. We have 
obtained minimal stability conditions which are necessary and sufficient in the Markovian 
case. Stability analysis also covers arbitrary initial state of the system. The conditions 
have clear probabilistic interpretation and can be easily applied to a number of important 
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particular cases. Examples of such particular cases are Erlang model and deterministic 
service model. We have observed that the stability condition for the Erlang model with 
retrial is insensitive to the distribution shape. 
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